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ABSTRACT 

We identify a new mechanism of internal dissipation of rotational kinetic 
energy in spinning dust grains, arising from the reorientation of nuclear angular 
momentum, e.g., spins of protons. Grain rotation induces magnetization of the 
nuclear spin system, with net alignment of nuclear spins parallel to the grain 
angular velocity. When the grain does not rotate around a principal axis, the 
nuclear magnetization vector precesses in grain body coordinates, resulting 
in dissipation of energy. The analogous process involving electron spins was 
discovered by Purcell and termed "Barnett relaxation" . We revisit the physics 
of the "Barnett relaxation" process and correct the estimate for the Barnett 
relaxation rate. We show that nuclear relaxation can be orders of magnitude 
more important than Barnett relaxation. This finding implies that the processes 
of "thermal flipping" and "thermal trapping" are important for a broad range 
of grain sizes. 



Subject headings: ISM: Atomic Processes, Dust, Polarization 



1. Introduction 

The polarization of starlight, first discovered by Hiltner (1949) and Hall (1949), 
implies that the axes of grains are partially aligned with the interstellar magnetic field. 
Alignment of a rotating interstellar grain, whatever its cause may be (see a list of alignment 
mechanisms in Lazarian, Goodman & Myers 1997), requires the alignment of the grain's 
body axes with its angular momentum J. Purcell (1979) noted that internal relaxation 
of energy within a rotating grain should cause alignment of the grain's axis of maximal 
moment of inertia with J. In the same pioneering study Purcell calculated the timescales 
for internal relaxation and found them to be much shorter than the gaseous damping time. 



1 This paper is dedicated to the memory of E.M. Purcell. 
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He discovered the process that he called "Barnett relaxation" and identified it as apparently 
the dominant mechanism of internal relaxation. 

Barnett relaxation is associated with the Barnett effect, i.e., the spontaneous 
magnetization of a rotating paramagnetic body (see Landau & Lifshitz 1962). A 
paramagnetic body has unpaired electron spins and the Barnett effect can be understood 
in terms of a rotating lattice sharing its angular momentum with the electron spin system. 
As a result, the rotating body becomes magnetized with magnetization M antiparallel to 
the angular velocity fl. For a freely rotating body fl and therefore M precess about J 
in inertial coordinates; in body coordinates fl and M precess around the body's axis of 
maximal moment of inertia (see Purcell 1979). This causes time- varying magnetization of 
the grain material and therefore entails dissipation. 

Spitzer & McGlynn (1979) related the Barnett relaxation time scale with the degree of 
disorientation of suprathermally (i.e. much faster than thermally) rotating grains during 
crossovers. The effect of suprathermal rotation was discovered by Purcell (1979), who 
identified recoils from H 2 formation events, occurring at catalytic sites on the grain surface, 
as the dominant mechanism of suprathermal spin-up.0 

The theory of crossovers was revised recently (Lazarian & Draine 1997, 1999; hereafter 
LD97, LD99) to include thermal fluctuations within the grain material (Lazarian 1994, 
Lazarian & Roberge 1997). The main conclusion by LD97 is that grains with radii 
a > a c ~ 1.5 x 1CT 5 cm experience only minor disalignment during crossovers if Barnett 
relaxation is the dominant dissipational mechanism. As grains are essentially immune from 
disorientation during their phases of suprathermal rotation, during which paramagnetic 
relaxation (Davis Sz Greenstein 1951) aligns them with the magnetic field, the modest 
disalignment during crossovers entailed nearly perfect alignment of a £ 1.5 x 10 _5 cm 
interstellar grains in diffuse clouds with the interstellar magnetic field. 

LD99 showed that for grains smaller than a c the crossovers are mostly performed via 
the process of "thermal flipping" and this would entail a high degree of grain alignment if it 
were not for the associated process of "thermal trapping"; the latter occurs when thermal 
flipping is so rapid as to effectively average out the systematic torques due to H 2 formation, 
etc. We note that the value of a c was calculated in LD97 and LD99 assuming Barnett 
relaxation to be the dominant mechanism of internal energy dissipation. In the presence 
of more effective relaxation mechanisms the coupling between rotational and vibrational 
degrees of freedom would be larger thermal flipping should be facilitated, and the value of 



2 A competing mechanism of spin- up due to radiation torques was recently proposed by Draine & 
Weingartner (1996, 1997). 
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a c would be increased. 

In this paper we show that the dynamics of the nuclear spinsQ leads to an important 
new process of internal dissipation in a rotating grain. Although the role of nuclear spins for 
paramagnetic relaxation was mentioned in the classic study by Jones & Spitzer (1968), and 
their role in developing a Barnett moment in a non-paramagnetic substance was pointed 
out by Purcell (1979), the internal relaxation associated with these spins was overlooked. 

Below we describe this new magneto-mechanical effect, which we term "nuclear 
relaxation" . We show that under a broad range of conditions this process dominates the 
dissipation of rotational kinetic energy in a rotating grain. As a result, "nuclear relaxation" 
plays an essential role in the dynamics and alignment of interstellar grains. 



2. Dynamics of Nuclear Spins 

Purcell (1978) noted that an analog of the Barnett effect exists for nuclear spins.0 
If a rotating body has initially an equal number of nuclear spins directed parallel and 
anti-parallel to the angular velocity f2, it can decrease its kinetic energy, at constant total 
angular momentum J, if some of the angular momentum is transfered to the nuclear spin 
system. Increasing the projection of the nuclear angular momentum along J by +K (at 
constant J) reduces the rotational kinetic energy by hfl. If the rotating body is allowed 
to come into thermal equilibrium (without exchanging angular momentum) with a heat 
reservoir of temperature Td us t then particles of spin S develop a net alignment per particle 

m exp(mhQ/ fcTdust) ^ 
EL-s exp(mhfl/kT dnst ) 

Note that this does not depend on the magnetic moment /i. 

In Table 1 we list a number of non-zero spin nuclei which could be important for 
interstellar grains. Carbonaceous grains will presumably contain 13 C at the ~ 1% level 
(n ~ 10 21 cm -3 ), and 1 H at perhaps the 10% level (n ~ 10 22 cm -3 ). Silicate grains 
will contain 29 Si at the ~ 1% level (n rs 10 21 cm -3 ), and could also contain comparable 
abundances of 57 Fe, 27 A1, and 55 Mn, as well as 1 H. For purposes of discussion we will assume 
the nuclear spin system to consist of protons with n n = 10 22 cm -3 . 



3 For brevity, we will use the term "nuclear spin" to refer to the total nuclear angular momentum Ih. 

4 This effect could be called the "nuclear Barnett effect" as opposed to the ordinary or "electron Barnett 
effect". 
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As the number of parallel and antiparallel spins becomes different the body develops 
magnetization. The relation between Q and the strength of the "Barnett-equivalent" 
magnetic field ii^E (Purcell 1979) that would cause the same nuclear magnetization (in a 
nonrotating body) is given by 

Hg = — n , (2) 

where g n is the so-called nuclear (/-factor (see Morrish 1980), and fijss = eh/2m p c is the 
nuclear magneton, smaller than the Bohr magneton by the electron to proton mass ratio, 
m e /m p . 

Consider a classical picture that would correspond to the concept of the Barnett- 
equivalent magnetic field. Imagine the (proton) spins as small positively charged tops 
within the rotating body. In body coordinates their motion is given by 

§| = *.xn . (3) 

where f2 is the angular velocity of the body in an inertial frame, and the prime signifies 
differentiation with respect to the rotating frame. It is evident that this equation is 
similar to the equation of motion in an inertial frame if a "Barnett-equivalent" magnetic 
field defined by eq. (0) were present. In other words, the motion of spins in the rotating 
body-coordinate frame is as though a fictitious "Barnett-equivalent" magnetic field 
were present in a nonrotating body. 



3. Nuclear Susceptibility 

The zero-frequency paramagnetic susceptibility of the nuclear spins with magnetic 
moments fi n = <? n A*N is given by Curie's Law (Morrish 1980) 

We wish to estimate the time scale on which the net alignment of the nuclear spin system 
will respond to changes in the direction of fl. When an oblate grain rotates freely about 
an axis that does not coinside with its principal axis of largest moment of inertia a, fl 
precesses about J. Due to the Barnett effect the magnetization of nuclear spins should 
follow fl. Interactions within the nuclear spin system (i.e., coupling of one nuclear spin to 
the magnetic field arising from other nuclear spins) conserve spin angular momentum, and 
therefore cannot change the direction of magnetization. However, the component of fl along 
J is constant, and so even very slow exchange of angular momentum between the lattice 
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and the nuclear spin system will result in magnetization of the nuclear spin system along 
J. Since J precesses in body coordinates, even magnetization along J will have dissipation 
associated with it. 

Consider first the coupling of the nuclear spins to the electron spin system.^ The latter 
is coupled to the lattice on the "spin-lattice" timescale r s i. If we assume that the electron 
spin system reacts relatively quickly to the change in rotational velocity J?, then the nuclear 
spin system will become aligned by exchange of angular momentum with the electron spin 
system. 

The "internal" magnetic field due to the electron spin system is Hi pa 3.8n e [i e (van 
Vleck 1937), where n e is the density of "unpaired" electrons, and /i e pa fi B is the magnetic 
moment per unpaired electron, where /zb = eh/2m e c is the Bohr magneton. Because the 
electrons themselves precess with characteristic frequency u e = g^^Hi/h, where g e ~ 2 
is the electron "g-factor", the correlation time for Hi is of order The fluctuating 

magnetic field causes the nuclear spins to undergo a random walk; each step lasts ~ w" 1 , 
and the change in direction of the nuclear moment per step is 59 ~ (g^^Hi/h)^' 1 <C 1. 
The initial phase relationship is lost after a time r ne pa (1/ '59) 2 uj~ 1 . Therefore the time for 
the nuclear spins to decorrelate with their initial direction - which is also the time scale on 
which they "relax" to the new alignment of eq. ([!]) - is 



The nuclear spin system can exchange angular momentum with the electron spin system on 
the timescale r ne . 

Curiously enough a quicker nuclear spin-spin relaxation arises from the direct 
interaction of nuclear moments. The magnetic field of the neighboring nuclei is ~ 3.8n n /i n , 
and the nuclear spin-spin relaxation time is 



For grains with high percentage of 1 H, r nn may be a factor of several smaller than r ne . The 
net nuclear relaxation rate r^ 1 = r m } + r~^. 

For an applied magnetic field of constant magnitude, rotating at frequency u, the 
magnetic susceptibility contributed by the nuclear spin system can be estimated to be (see 



5 We use the term "electron spin system" to refer to the system of paramagnetic ions and electrons. 
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Draine & Lazarian 1999) 



XnH « Xn(0)- : (7) 

(1 — lUT n /2y 



X >) = Im[ Xn M]= Xn (0)- ^" 



[1 + (ut d /2] 



212 



4. BARNETT RELAXATION REVISITED 



Purcell (1979) assumed that the dissipation arising from the Barnett effect in a rotating 
grain is the same as would occur if the Barnett equivalent magnetic field were rotating in a 
stationary grain. The analogy, however, is not complete. Consider a freely rotating oblate 
grain with eigenvalues of the moment of inertia tensor In, I±, I±, with I\\ > I±. Let a be the 
principal axis of largest moment of inertia, In. If the grain is not rotating about a principal 
axis, then in inertial coordinates, a and the angular velocity fl each precess around the In 
grain body coordinates, the angular velocity fl precesses around a; the frequency of this 
precession is 

/, „ (h — l)Jcos9j . . 

u 1 = (h- l)n cos 9q = ± L (9) 

(see Purcell 1979), where h = I\\/I±, Qq, is the angle between fl and a, and 6j is the angle 
between J and a.f] 

Let T e i ~ 10 _6 s be the electron spin-lattice coupling time. In the limit u>it c \ <C 1, the 
electron spin system will respond to the instantaneous angular velocity, with a magnetization 
X-^bej "where = hQ,/g e /j,^ is the Barnett-equivalent field for the electrons. In body 
coordinates, this magnetization has a component X-^be s i n which is rotating around a, 
and as a result there will be energy dissipation at a rate dE/dt = ^{oJi)V {H^) 2 sin 2 9q. 

In the opposite limit 0JiT e \ ^> 1, the electron spin system is almost decoupled 
from the lattice. The electron spins interact magnetically, but the total spin of the 
electrons is conserved, and cannot follow fl as it precesses around J. However, since the 
spin lattice coupling will not be zero, eventually the electron spin system will become 
magnetized in response to the (stationary) component of fl along J, with a magnetization 
X-Hbe cos (8n — Oj)- This magnetization has a component X-^BE cos ($n — #/) sin 6*j 
perpendicular to a, and in body coordinates this component of the magnetization will 
rotate at frequency U\. We can write the energy dissipation rate as 

dE/dt = x'^ViH^l) 2 sin 2 ^/(a;!, 9j) (10) 



6 The two are related via tanfc = ht&nOj. 
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where 



/ 



1 + (c^Td) 2 sin 2 9j cos 2 (^ - 0j)/ sin 2 9 Q 



1 + (^ir cl ) 2 

gives the asymptotic behavior for u>iT e \ <C 1 and uj\T e \ 1. 

The kinetic energy of the grain is (J 2 /2/||)[(l + (h — 1) sin 2 6*j]. Thus, 
dE/d6j = Ju 1 sin9j, and Barnett relaxation entails 



in) 



^ = -o-i>^, • , 

at lull \g n ^J \vi 



h ' ( ^sm0j<xya0jf{h,0j), 



(12) 



The Barnett relaxation is slower than estimated by Purcell (1979) by the factor / < 1. 
However, this factor is of order unity, and we will take / = 1 in the remaining discussion. 

For an ax a^/3xa^/3 brick, Barnett relaxation gives 



r BR ^2.1 x 10 8 pV 5 



(13) 



a = a 



/10 5 cm and Jd = [I\\kT d /(h — l)] 1 ^ 2 is a characteristic 



where p = p/(2gcm 
angular momentum appearing in the theory of crossovers (LD97, LD99). 



10 1 



10 1 



"?io 8 



10 6 r 



10 4 f 




1 10 r____ — — 



o 10 B r 
10 6 r 



crossover timet c y_ — _ — — 



(J/J d ) 2 =l 



Nuclear Relaxation 



10 4 r 

io- 6 



Barnett Relaxation 

_J 1 1 1 1 I I I I 1 L 



_J 1 1 1 | | 



io- 



a(cm) 



io- 4 



Fig. 1. — Time r for alignment of a with J for a freely rotating grain, due to dissipation 
associated with nuclear spins (nuclear relaxation) and electron spins (Barnett effect). Results 
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are shown for grains with (J/J d ) 2 = 1 (lower panel) and 10 (upper panel). Also shown is 
the "crossover time" t c = J/ J», where the torque J» is due to H 2 formation with a density of 
active sites 10 13 cm -2 cm -2 . The results shown here are for an axbxb "brick", with b/a = \/3, 
spinning with cos6j = 0.5. The grain is assumed to have p = 2 gem -3 , unpaired electrons 
with n e = 10 22 cm -3 , and protons with n n = 10 22 cm -3 . 



5. Nuclear Relaxation In a Rotating Grain 



Nuclear relaxation is more effective than Barnett relaxation. The time for nuclear 
relaxation to align a with J is 

2 



T"NR = 



T T 2 n 2 n 2 
1 a 1 ± i/nPN 



(h-l)VJ 2 k 2 Xn(0)r n 



1 + 



U)iT n 



For an axbxb brick with b/a = \^3, we have h — 1 — 1/2, and 



Jd\ fn e \ fg n \ 4 [2.7(1 



2 r 



fin 



1 + 



1 2 



(14) 



(15) 



Figure 1 shows Tnr, Tbr and the overall relaxation time r = (r^ + Tg^) -1 as a function 
of grain size, for grains rotating with J/ J d = 1 and J = lO 1 / 2 ^. For slow rotation rates 
(large a) the timescale tnr is ~ 10 6 times shorter than tbr! As a result the critical grain 
size for which thermal fluctuations dominate the dynamics of crossovers is increased by an 
order of magnitude. We conclude that "thermal flipping" (LD99) must therefore be an 
essential part of crossover dynamics for all grains with a ^ 10~ 4 cm. 



6. Discussion 

The first striking question is why an effect so feeble as nuclear magnetism can be so 
important in terms of internal relaxation. The answer is that spins align in a rotating 
body because of their angular momentum, not because of their magnetic moments. The 
magnetic moment enters only as a means for the spins to exchange angular momentum 
with the lattice. The coupling within the electron spin system is very effective, with the 
result that there is minimal "lag" in the electron spin alignment when the grain angular 
velocity S2 changes in grain body coordinates. In the case of the nuclear spin system, 
however, the value of the nuclear magnetic moment is large enough to provide significant 
spin-spin coupling, but weak enough so that there is a significant lag in the nuclear spin 
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alignment when 17 precesses around the grain axis a: the coupling is "just right" for nuclear 
relaxation to be extremely effective for 10~ 5 £ a £ 10~ 4 cm grains. The efficiency of nuclear 
relaxation drops for sufficiently high frequencies and therefore Barnett relaxation dominates 
for a 5 x 10~ 6 cm, as seen in Figure [T]. 

LD99 found that "thermal flipping" was an important element of grain dynamics for 
sufficiently small grains. Assuming the Barnett effect to dominate internal dissipation, 
they estimated that thermal flipping would be important for a £ 10~ 5 cm grains, and they 
speculated that for a £ 5 x 10~ 6 cm grains, thermal flipping would be so rapid as to interfere 
with the ability of systematic torques to spin these grains up to suprathermal rotation rates. 

In the present paper we have shown that internal relaxation associated with nuclear 
spin alignment can be many orders of magnitude more rapid than due to the Barnett 
effect for a £ 10~ 5 cm grains. As a result, the phenomena of thermal flipping and thermal 
trapping become important for grains as large as ~ 10 -4 cm. 

7. Conclusions 

The results of this study can be summarized as follows: 

1. Freely rotating interstellar grains are subject to "nuclear relaxation" when their 
angular momentum J is not along a principal axis. 

2. For slowly rotating interstellar grains the dissipation due to nuclear relaxation is 
many orders of magnitude higher than that due to Barnett relaxation. 

3. The critical size a c below which thermal fluctuations are effective during a crossover 
increases to ~10~ 4 cm. Thermal flipping and thermal trapping must therefore be essential 
elements of the dynamics of all a £ 10~ 4 cm interstellar grains. 

We thank R.H. Lupton for availability of SM. This research was supported in part by 
NASA grant NAG5-7030 and NSF grant AST-9619429. 
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Table 1. Nuclear Properties 3, 



Nucleus 


J b 


Abund.(%) 


/i n //iAT C 


9n 




1/2 


99.998 


2.675 


3.089 


13 C 


1/2 


1.108 


0.673 


0.777 


14 N 


1 


99.635 


0.193 


0.137 


17 


5/2 


0.037 


0.363 


-0.123 


35Q 


3/2 


75.53 


0.262 


0.135 


27 Al 


5/2 


100 


0.697 


0.236 


31p 


1/2 


100 


1.083 


1.250 


57p e 


1/2 


2.19 


0.086 


0.100 


29 Si 


1/2 


4.70 


0.531 


-0.614 


55 Mn 


5/2 


100 


0.660 


0.223 



a From Robinson (1991) 

^total angular momentum quantum number 

c /iN = eh/2m p c — 5.05 x 10~ 24 erg/gauss. 



